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Comparison of Laser Models 

Notation. 

In this report the following notations will be used: 


%\ = [Ri] 
x 2 = [R] 
x z = [Ra] 


*4 = [/a] 

*5 = [/*] 

X <5 = [I\ 
17 = [p\ 


where [ ] denotes concentration and R is a perflouride radical, I is iodine, /* is an excited 
state of iodine, and p is the photon density. 

Oscillatory Model. 

The following is the unsealed version of the time varying chemical kinetics equations 
for the simulation of the iodine laser operation 


dx i 

IT 
dx<2 
IT 

- kr? 

dt 

d%4 

If 

- = -kiX 2 X& - cixixsxe - C3X4X5X6 - Qixixt - <? 2 x 4 x 6 


= fciZ 2 X 5 + *2*2 *6 ~ tMl - *4*1*2 
= -A1X2Z5 - *2*2*8 - 2^3X5 + ^1*1 _ ^4X1X2 
= * 3*2 + A4X1 X2 

= CiXiX 5 Xfl + C2X1XJ -1- C3X4X6X6 - if) 2X4 + C4X4X3 


- Az 5 + ^1X1 + 0.51^2*4 - T, 


max 


d%$ 

= 1.49^2X4 + Q1X1X5 + Q 2 *4*6 + Az 5 - CiXiXgXe + T max 


1 * 4 X j *2 ^2^1 C3X4X3X3 — * *3X3X3 2C4X4ZQ 


dx 7 

IT 


\jrrr (l c ) 


where I'm** = 


n , log(R 1 R 3 ) 

A max i 

cx 7 (x 8 - ix«) 


„ /0.18V 

X7 + 2 (-ET j 


AX 6 


Aoo + Roo*i 

The above equations are solved subject to the initial conditions 


Xi (0) = P * (3.5R16) and x,(0)=0, t = 2,...,7 
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where P is the pressure in torr. The above system was solved using a variable step size 
Runge-Kutta-Fehlberg alogorithm. The figures 1.1 through 1.4 illustrated a typical set of 
solution curves to the above system of differential equations. Note the violent oscillations 
which occur during start-up of the laser. These oscillations require that small step sizes 
be taken during the numerical integration process in order to achieve some desired degree 
of accurracy. The numerical method becomes unstable and the results are meaningless if 
too large a step size is selected. This is not an uncommon situation in trying to numerical 
solve chemical kinetic systems of equations. 

Nonoacillatory Model. 

The inversion term (/*] - ^ (/] = £5 - 5X6 which occurs in the term as well as 
the equations which define the rate of change of [/*] = Z5 and [/] = xe, is the term which 
causes the oscillations in the nonlinear system of coupled ordinary differential equations. 
Under quasi-steady state conditions we set the rate of change of this term equal to zero. 
That is we require that 

3 (i j '> - H = 0 

and then solve for r max which is then relabled as (T mnx ) ta . We find that 

(r m o*) = — Ax$ — Q2X4X5 + {2^1X1 — 0.47^2X4 — 2A1X2X5 + £3X3X9 
- C1X1X5X6 + 2 C 3 XixJ - C3X4X5X6 + 2 c 4 X 4 X^ - k 4 x iX 2 }/3 

Initially the r max term is removed from the above equations as it is assumed to be initially 
zero. As time increases, (/*] and [i] build up and when the threshold condition is met we 
have , . _ _ , 

Here it is assumed that steady state conditions have been achieved and gains balance with 
losses. The above steady state value is subsituted for r ma *. In this way the ocillatory 
character of the solution is removed and the integration can be made to proceed at a much 
faster rate. The results from this second nonoscillatory model were generated by L.Stock 
and are presented in the figures 1 . 5 , 1 . 6 . Observe how the oscillatory character of the 
solutions has been removed. Also, as a check on the model, the values of the power output 
were compared for both models and both values were found to be in agreement. The figure 
1.7 lists the values of the reaction rate coefficients used in both models. 

Currently there are four models for the simulation of an iodine laser. These models 
can be described by: 


1. A time dependent model which is given by the above set of coupled nonlinear differ- 
ential equations which describe the chemical kinetics at a fixed point as a function of 
time t. 

2. A quasi-steady state model whereby the r max term is replaced by its steady state 
value. This model is also a time dependent model with much faster integration times. 

3. A noncompressible model which solves the above system of differential equations as 
a function of distance z along the axis of the laser tube. This model is characterized 
by substituting a material derivative in place of the time derivatives in the above 
differential equations. For an assumed steady flow, we make the substitutions: 

in = 211+ 211* 

dV dt dz dt 

and assume that ^ = uj is a constant. The time dependence is removed by assuming 
steady state values. Problems of stiffness arise in the numerical integration of this 
system. 

4. A compressible flow laser model which solves the above system of differential equa- 
tions with the material derivative replacing the time derivatives. In addition to these 
equations there is added: (i) a continuity equation (ii) an equation of state (iii) a 
momentum equation and (iv) an energy equation. 

The figures 1.8 illustrates power vs beam radius which was obtained from the model 4, 
above assuming a parabolic velocity profile. The figure 1.9 illustrated power vs reflectivity 
for various pumping distances along a 45cm laser tube. 

Current research is being directed toward: 

1. Parameter studies using the compressible flow laser model. 

2. Development of a two pass amplifier model. 

3. Solve system of equations describing operation of high powerhh iodine MO PA (Master 

Oscillator Power Amplifier). These equations are derived from the Frantz-Nodvik 

Amplifier theory. 
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Fig. 1.3 Inversion Term vs. Time 
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Fig. 1.4 Power vs. Time 
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Power vs. Time 


TABLE L - Reaction rata coefficients and their associated reactions for 
the lasant n-Cjy along with the parameters used to modify the optical 
time constant (ref. 2). 



Lasant 

n~C 4 F 7 I 

Reaction Rate Coefficient 


Reaction* 

Symbol 

(cm 3 )"/*® 6 

R + 

r ♦ ri 

K, 

.5600 X 10"“ 

R + 

1 - RI 

k 2 

.2300 X IQ " 10 

R + 

R - R t 

K, 

.2600 X 10 ' 11 

R + 

RI - R t + 1 

K* 

.3000 X 10 " 14 

R + 

l 2 - RI + 1 

K, 

.0000 X 10** 

1 * + 

1 + RI - l a + RI 

c, 

.3200 X 10-“ 

1 + 

1 + RI ■» l 2 + RI 

C. 

.6500 X IQ " 81 

1 * ♦ 

1 + 1 , - l 4 + 1 , 


.6000 X 10 -51 

1 ♦ 

1 * l 2 * l 2 + I 2 

c* 

.3600 X 10-* 

r + 

RI - 1 * RI 

0 ,' 

.2000 X 10 " 14 

i* + 


0 2 

.1900 X lO * 10 



T s 

1 .0000 



a 

.0000 


Fig. 1.7 Table I. Reaction Rate Coefficients 
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Fig. 1.8 Power vs. Beam Radius 



45 CM LASER WITH VARIABLE PUMPING DISTANCE AND VARIABLE REFLECTIVITY 
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